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We calculate, within the pseudopotential approximation, 
a one-dimensional scattering amplitude and effective one- 
dimensional interaction potential for atoms confined transver- 
sally by an atom waveguide or highly elongated "cigar" - 
shaped atomic trap. We show that in the low-energy scat- 
tering regime, the scattering process degenerates to a total 
reflection suggesting an experimental realization of a famous 
model in theoretical physics - a one-dimensional gas of im- 
penetrable bosons ("Tonks" gas). We give an estimate for 
suitable experimental parameters for alkali atoms confined in 
waveguides. 



Rapid progress in producing Bose condensates of al- 
kali atoms Q opened up new areas of ultra-low energy 
collisional physics. Concurrently, work has progressed 
on confinement of atoms in the light-induced and 
magnetic-field-induced |3| atom waveguides. 

We develop a theory for the binary atomic collisions in 
presence of a transverse external confinement. Using the 
pseudopotential approximation, we derive an expression 
for an effective one-dimensional scattering amplitude and 
show that the interparticle interaction can be approxi- 
mated by an effective one-dimensional 5-potential of a 
known strength. In the case of a dilute atomic gas, when 
the three body collisions are negligible, our effective po- 
tential can also be used to describe quasi-one-dimensional 
many-body systems: a project on the experimental real- 
ization of the one-dimensional Bose condensate is already 
presented in the literature . Results of our paper will 
allow one to properly take into account the trap-induced 
corrections to the strength of the atomic mean-field po- 
tential. 

Furthermore, we show that in the low-energy scatter- 
ing limit (fc z |aiD| *C 1), the effective one-dimensional 
scattering degenerates to a total reflection (hk z is the lon- 
gitudinal component of the atomic momentum and Oid 
is the one-dimensional scattering length defined below). 
This conclusion allows us to suggest an experimental re- 
alization of another famous model in theoretical physics 
- a one-dimensional gas of impenetrable bosons M also 
referred to as a "Tonks" gas. Such a system provides 
us with an unusual example of a boson-fcrmion dual- 



ity: the elementary excitations of such a bosonic sys- 
tem obey Fermi statistics || (which is possible only in 
one dimension 0). The Tonks gas is, therefore, a sys- 
tem complementary to the one-dimensional Bose conden- 
sate: the latter requires the high-energy scattering regime 
(&z|oid| ^ 1) H an d its excitation spectrum is repre- 
sented by the Bogoliubov bosons. 

To describe the binary collisions between cold atoms 
confined in a waveguide, we suggest the following model: 

(a) The waveguide potential is replaced by an axially 

symmetric 2D harmonic potential of a frequency 
ui±. The forces created by the potential act along 
the X-Y plane; 

(b) Atomic motion along the Z axis is free; 

(c) Interaction between the atoms is modeled by the 

Huang's pseudopotential 



U(r) = gS(r)^(r-) , 



(1) 



where g = 2Trh 2 a/fi is the potential strength, a is 
the s-wave scattering length for the "true" interac- 
tion potential, ^ — m/2 is the reduced mass, m is 
the mass of the atoms. The regularization operator 
■ftp (r • ) , that removes the 1 jr divergence from the 
scattered wave (see |!|), plays an important role in 
the derivation below; 

(d) Atomic motion (both transverse and longitudinal 
components) is "cooled down" below the transverse 
vibrational energy ?iuj±. We will justify this condi- 
tion in the second subsequent paragraph. 

The harmonic nature of the confining potential allows 
the separation of the center-of-mass and relative motions. 
The Schrodinger equation governing the relative motion 
reads 



~2 o 

Ez-+g6(r)— (r-) + H ± (p x ,p y ,x,y)}* = EH> . (2) 



where r = r2 ri is a relative coordinate for atoms 1 
and 2 and 



Px+Py pLo\{x 2 +y 2 ) 
2[i 2 



(3) 
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is a 2D harmonic oscillator Hamiltonian. 

We will suppose that (i) the incident wave 
e lk * z <f)n=o,mz=o(p) corresponds to a particle in the ground 
state of the transverse Hamiltonian @ (p = ^ x 2 + y 2 ) , 
and (ii) the longitudinal kinetic energy of the incident 
wave is limited by the energy spacing between the ground 
and first axially symmetric excited state: 
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(4) 



Here B„ imi = htu±(n + 1) is the energy spectrum of the 
2D harmonic oscillator, n = 0, 1, 2, . . . , oo is the princi- 
pal quantum number, m z is the angular momentum with 
respect to the Z axis, m z = 0, 2, 4, . . . , n (1, 3, 5, . . . , n) if 
n is even(odd). The asymptotic form of the scattering 
wave function then reads 
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+ /oddSign(z)e ifc -l z l } <£ ,o(p) , (5) 



where the first term in the curly brackets represents the 
incident wave, the second and third terms give the even 
and odd scattered waves respectively, and f eV en(k z ) and 
/odd(fcz) are the one- dimensional scattering amplitudes 
for the even and odd partial waves respectively. Note that 
the transverse state (n = 0, m z = 0) remains unchanged 
after the collision. Transitions to the higher (n > 0) 
modes are forbidden due to either angular momentum 
or energy conservation (if the condition (Q) is satisfied). 
We would like to stress however, that during the collision 
itself a virtual excitation of the high energy axially sym- 
metric modes (n > 0, m z = 0) is not forbidden and it is 
properly taken into account below. One of the goals of 
this paper is to calculate the amplitudes (^). 

For the zero-range potential ([!]) the one-dimensional 
scattering amplitudes (|5|) can be calculated analytically. 
To perform this calculation we expand the wave func- 
tion ^(z,p) to a series over the eigenstates of the trans- 
verse Hamiltonian (^), substitute the expansion into the 
Schrodinger equation (Q) using E = h 2 k 2 z /2p, + huj ± for 
the energy in the right hand side, and then apply the 
asymptotic conditions (^|) along with the conditions of 
the continuity of the wave function and its derivative. We 
obtain the following expression for the scattering ampli- 
tudes: 



fe 
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where rj is the r — ► limit of the regular (free of the 1 jr 
divergence) part of the solution \E': 
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(J) 



The expression for the wave function reads 



*(2,p=0) 
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exp(ifc 2 z) 
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(8) 



where the function A is defined as A[£, e] = 
J2^>=i ex P( — V s ' + e £)/V s ' + 6 (ll); the sum over s' — 
n/2 originates from a sum X^^Lo 2 4 ••■ over the trans- 
verse states n ,m^=o, and a± = {%/ pwjj 1 / 2 is the size of 
the ground state (j)n=o,rn z =o of transverse Hamiltonian 
(^). Above we have used a simple relation \4> n , mz =o(p = 
0)| 2 = \j-Ka\ for the wave functions of the 2D harmonic 
oscillator. The value </>o,o(p = 0) has been chosen to be 
real and positive. 

The regular part r\ of the wave function is not explicitly 
defined yet, being a solution of the equation ([?]). Note 
that the order of the partial derivative d/dz and the 
sum over s' in the A-function can not be interchanged 
because the series A does not converge uniformly as 
£ — > 0. The value of r\ can be extracted from the equa- 
tion (0) using the following expansion for the function A: 
A[£, e] = 2ji + £(e) + £i(e)£ + . . . , where the zero-order 
term of the expansion has a form £(e) = — C + £(e) and 



(9) 



/ -= - V -= = 1.4603... 

Jo V^ 7 ftx v^y 

n=l 

To prove this formula one should subtract and add a sum 

E"=i Is>-1 ds " exp(-V^ 7 £)/\ / s 77 = 2/£ to the function 
A. The sum of the differences between the A terms and 
the terms of the above sum converges uniformly and the 
derivative in (Q) can be easily calculated. Here t[£] is the 
Riemann zeta-function. 

We now write down the final expression for the one- 
dimensional scattering amplitude (||): 



/even(^) 



1 + ik z a 1D - (ik z ax/2) £(-fc>i/4) 



(11) 



o((fe«i) 3 ) 



where the function £(e) is given by (|io|). Here 
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(12) 



is the one- dimensional scattering length and the constant 
C is given by (^|). In analogy with three-dimensional 
scattering, the one-dimensional scattering length is de- 
fined as a derivative — dA/dk z \k z ^o+ of the even- wave 
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scattering phase. The scattering phase A(fc z ) is defined 
through the even solution 'F oc sin[fc z |z| + A(k z )](po,o{p) 
of the Schrodinger equation (||) . The formula for the one- 
dimensional scattering amplitude ( |ll| ) is the key result of 
this paper. 

The expression ( pi] ) is valid for any strength of the 
transverse confinement. Note however, that the tight 
confinement limit a± <C \a\ makes sense only if the s- 
wave scattering amplitude for the "true" 3D finite-range 
interatomic potential (approximated in this paper by the 
Huang's pseudopotential (|IJ)) shows a zero-energy res- 
onance; namely if the s-wave scattering length a (the 
same for both potentials) is much greater than the effec- 
tive range tq of the "true" potential JhJ . Indeed the 3D 
pseudopotential approximation (jl|) is only valid for ve- 
locities which are lower then the inverted effective range 
(k -c r^ 1 ) (see E2]), and therefore requires a lower bound 
for the transverse size: aj_ ~ kj 1 tq. This condition 
and the tight confinement regime criterion (aj_ <C \a\) 
are consistent only in the resonant case: \a\ 3> tq. 
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FIG. 1. Transmission probability as a function of the in- 
cident momentum. Solid lines correspond to the exact ID 
scattering amplitude ([H]) . Dashed lines correspond to the ID 
(^-potential approximation /even ~ /oven- For a±/a = ±10 the 
result of the 5-potential approximation almost coincides with 
the exact amplitude and is not shown. 



For low velocities the exact scattering amplitude 
( |ll| ) can be approximated by a scattering amplitude 
/even(^z) = — 1/(1 + ik z aij)) for a one- dimensional 8- 
potential 




Um{z) = gir>S(z) 
of a coupling strength 

Sid = - = 9 |0o,o(O)| 2 (l - C — 

A* a iD \ a_L 



(13) 



(14) 



where |0o,o(O)| 2 = \ji^a\. To illustrate this state- 
ment we plot (FigJ^) the transmission coefficient T = 
|1 + /oven + /odd | 2 calculated using the exact scatter- 
ing amplitude (|ll|) along with the results of the ID 5- 
potential approximation f cvcn s» /,f vcn . Recall that in 
both cases the odd- wave scattering amplitude vanishes: 

/odd = / dd = 0- 

The effective potential ( |l3| , |T4[ ) can be shown to repro- 
duce the low-energy scattering properties of the radius 
a hard spheres in presence of a transverse trap. Recall 
that the Huang's potential ([!]) plays exactly the same 
role in the free space. By analogy with the "free-space" 
case ||, the the limits of applicability of our potential 
can be extended to the many-body problems to describe 
interactions between the atoms confined to the ground 
transverse state. 

Note that the one-dimensional gas of bosons interact- 
ing via a ti-potential is already widely investigated in the- 
oretical physics fijl . One of the most interesting models 
belonging to this class is a one- dimensional gas of im- 
penetrable bosons. Formally this model corresponds to 
the infinitely strong repulsive interaction between atoms: 
.9id - * +oo. More rigorously the impenetrable bosons 
regime corresponds to the low-energy scattering limit 



k z \a 1D \ < 1 



(15) 



when the corresponding transmission coefficient T ap- 
proaches zero (see Fig|l]). (Recall that for positive inter- 
action strength guy the corresponding one-dimensional 
scattering length aiD is negative.) A remarkable feature 
of the gas of impenetrable bosons is a possibility of a one- 
to-one mapping between the bosonic system and a gas of 
noninteracting fcrmions. 

For N impenetrable bosons confined in a periodic- 
boundary-conditions one-dimensional box of a length L 
the ground state of the system ^ b is given by an absolute 
value of the ground state of the iV-particle ideal Fermi 



gas: 



where 



lip/ 1 



H f (zi,z 2 , ...,zn) 



1 



det(e 
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kj — ~~T~j ^ t ^Fcrmii +&Feri: 



(16) 

(17) 
(18) 



and fepormi = n(N — 1)/L is the one-dimensional Fermi 
radius. In Fig.|| the zero-temperature one-body mo- 
mentum distribution w(k z ) for a system of impenetra- 
ble bosons in the thermodynamic limit is shown. (It is 
normalized as f^°° (dk z /kF etm i) w(k z ) = 1.) The dis- 
tribution has been calculated using the short-range and 
long-range expansions for the one-body spatial correla- 
tion function given in G3j. The shape of the peak at 
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the origin is given by w(k z ) w /Ooo\/fcFermi/27rfc z where 
Poo = ne 1 / 2 2- 1 / 3 A- 6 = 0.92418 . . ., and A = 1.2824 . . . 
is the Glaisher's constant. Note that for finite values 
of N and L the thermodynamic limit fails at small mo- 
menta k z ~ 2tt/L ~ 2k-p eTlrL i/N and the momentum den- 
sity curve is not valid in this domain. For comparison, 
we plot also the corresponding momentum distribution 
for an ideal Fermi gas. 

The system of impenetrable bosons may be realized in 
the atom wave guides Q combined with a longitudinal 
confinement between two potential barriers. The impen- 
etrable bosons regime may be identified by a presence of 
the 1 / \fk~ z peak in the momentum distribution of atoms 
(see Fig.f). Note that the low energy scattering condi- 
tion ( p.5[) puts an upper bound to the number of atoms. 
Indeed, using the fact that the maximal momentum of 
the relative motion between a pair of atoms is given by 
kmax — (p/m) max. j j/ (kj — kj>) — fcpermi one could show 
that the low energy scattering condition & m a£c |oid| <C 1 
leads to a limit for the number of atoms 



N -C N* = 



L 



7T OlD 



(X L LO±_ 



(19) 



where oid is given by Jl2]). For oj± = 2ir x 10 4 Hz and 
L = 3 cm the upper bound for the number of atoms is 
given by: N^ h = 5 x 10 3 for rubidium 87 (F = 2, mj? = 
+2, a = +110a Bo hr H|l); and N£ a = 6 x 10 2 for sodium 
(F = l,m F = -l,a = +52a Bo hr Q). Note that the 
above numbers can be improved in stronger traps since 
the bound N* increases approximately linearly with the 
confinement frequency. 
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FIG. 2. Momentum distribution for a system of impen- 
etrable bosons at zero temperature in the thermodynamic 
limit. Corresponding distribution for an ideal Fermi gas is 
shown for comparison. 



In conclusion, we calculate the one-dimensional scat- 
tering amplitude and effective one-dimensional interac- 
tion potential for atoms transversally confined by a two- 
dimensional harmonic potential. We suggest a realization 
of the Tonks gas - a one-dimensional gas of impenetrable 
bosons. From the experimental point of view the impen- 
etrable bosons regime corresponds to highly elongated 
traps (waveguides) at low temperatures (k^T <C htu±) 
and at low linear densities ( pa <C (2ir)~ 1 (a/a±) 2 , a > 0). 
We give an estimate for suitable experimental parameters 
for alkali atoms confined in waveguides. 
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